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Abstract 

We study diffusion in a network which is governed by non-autono- 
mous Kirchhoff conditions at the vertices of the graph. Also the dif- 
fusion coefficients may depend on time. We prove at first a result 
on existence and uniqueness using form methods. Our main results 
concern the long-term behavior of the solution. In the case when the 
conductivity and the diffusion coefficients match (so that mass is con- 
served) we show that the solution converges exponentially fast to an 
equilibrium. We also show convergence to a special solution in some 
other cases. 
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1 Introduction 

The aim of this paper is to study a non-autonomous dynamical system in a 
network subject to non-autonomous Kirchhoff knot-conditions at the vertices 
(knots). To be more precise we consider a finite, simple, connected graph 
G with m edges ei, . . . , e m and n vertices vi, . . . , v n . We study a diffusion 
system of the form 

^(t,x) = Cj (t)^(t,x) + Fj(t,x), (1.1) 

j = 1, . . . , m, t > 0, x € (0, 1). Here we think of Uj(t, x) as a function defined 
on the j-th edge e 3 - and impose that, for each t > 0, (ui(t, •),... ,u m (t, •)) 
is continuous on the graph, i.e. that the values of Uj(t, •) are compatible at 
the vertices. We further impose Kirchhoff conditions at the vertices 

m o 

E^(0^(i,v,) = 0. (1.2) 
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Here <1> = is the incidence matrix of the graph and Hj(t) are given 

positive conductivity factors. Finally we impose an initial value condition 

it(0,x) = u°(s). (1.3) 

We will show well-posedness and study the asymptotic behavior of the so- 
lution as time tends to infinity. 

Such evolutionary systems on networks have been studied by many au- 
thors, in particular in the autonomous case we refer to the monograph 
[LLS94] and the proceedings [ABN01]. Concerning the non-autonomous case 
we mention the work by von Below [Be88]. Here we use a different strategy 
than von Below, namely we use form methods to prove well-posedness. In 
the autonomous case this form method allowed in [KMS07] to establish a 
holomorphic semigroup, and the asymptotic behavior of the solutions was 
studied with the help of spectral theory. There are results on existence and 
uniqueness for non-autonomous forms, much developed by J. L. Lions, but 
the problem is to obtain solutions with sufficient regularity to identify the 
knot-conditions at the vertices. This is indeed possible with the help of a 
recent result [ADL012] which we use here. We need to assume that the 
conductivity factors )J.j(t) are Lipschitz continuous in time. But then we 
obtain a unique solution for each initial value u° G V which is continuous 
from R + with values in V where 

V = {/ € H l (0, l) m : / is continuous on the graph}. 

For studying the long-term behavior of the solution we need at first to 
show that the eigenvalues depend continuously on time. This is done with 
the help of a criterion which is known to specialists but is possibly nowhere 
formulated explicitly. We give precise information in the Appendix. 

This paper seems to be among the first studies of asymptotic behavior 
for non-autonomous evolution equation in a network. A recent manuscript 
[BDK13] treats a first order transport equation on the edges of a network 
with time- varying transmission conditions at the vertices using the approach 
of difference equations and evolution families. However, the nature of our 
problem and the techniques we use are completely different. 

It is natural that the asymptotic behavior depends on conditions on the 
coefficients. We consider several cases. If the diffusion coefficients Cj(t) and 
the conductivity coefficients fjij(t) are equal, then mass is preserved. We 
show in this case that the solution converges exponentially fast to an equi- 
librium (Section 5.1). As a second case we assume that bj(t) := satisfies 
a monotonicity condition (Section 5.2). Then the solution converges expo- 
nentially fast to a special solution. Finally, if bj(t) < cbj(t) (Section 5.3), 
then we can still prove exponential convergence, but we have less information 
on the limit solution. 
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The paper is organized as follows. Section 2 is devoted to abstract non- 
autonomous forms. We give the well-posedness result (Theorem 2.4) which 
we use later and prove continuity of eigenvalues as functions of t. In Section 
3 we explicitly describe our network diffusion problem and establish the well- 
posedness. In Section 4 we investigate when the solutions are positive. The 
qualitative behavior of solutions for t — > oo is studied in Section 5. Finally, 
we show in the Appendix that the A;-th eigenvalue of a self-adjoint operator 
is continuous with respect to convergence in the resolvent sense. 

2 Forms and associated operators 

In this preliminary section we introduce the abstract framework. Let V 
and H be two separable real Hilbert spaces such that V H, i.e. V is 

d 

continuously and densely embedded in H. We denote the scalar products 
and norms of H and V by (•,•)#> W'Wh an< ^ (v)y> ll'llv respectively. Let 
o: V x V — t- R be a bilinear and continuous mapping, i.e. 

\a(u,v)\ < M\\u\\ v \\v\\ v (u,v&V) (2.1) 

for some constant M > 0. We assume that a is H -elliptic, i.e. there exist 
constants a > and u> G R such that 

Rea(u, u) +uj\\u\\ 2 h > a||w||y (uEV). (2.2) 

If u = 0, we say that the form o is coercive. The operator A £ C(V,V') 
associated with the form a is defined by 

(Au, v) = a(u, v) (u,vEV). 

Here V denotes the dual space of V and (•, •) denotes the duality between 
V' and V. As usual, we identify H with a dense subspace of V' (associating 
every / € H with the linear form v \-¥ (f,v)u). Then V is a Hilbert space 
for a suitable scalar product. 

Seen as an unbounded operator on V' with domain D{A) = V the op- 
erator —A generates a holomorphic semigroup on V' . The semigroup is 
bounded on a sector if uj = 0, in which case A is an isomorphism. Denote 
by A the part of A on H, i.e. 

D(A) := {u G V : Au G H}, 
Au := Au. 

Note that in this case we have 

D(A) = {u G V : 3h G H s.t. a(u,v) = (h,v) H for all v G V}, 
Au = h. 
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Then the operator —A generates a holomorphic Co-semigroup on H (the 
restriction of the semigroup generated by —A to H). We call A the operator 
associated with a on H. Note that by the Lax-Milgram Theorem (A + 
A) : D(A) -> H is bijective and (A + A)" 1 G C{H) for all A > u. For more 
information, see e.g. the monographs [Ouh05, Chap. 1] or [Tan 79, Chap. 2]. 

We now consider a symmetric form a, i.e. we suppose that a(u, v) = 
a(v, u) for all u, v G V. The following estimates follow from the -ff-eilipticity 
of a. 

Proposition 2.1. Assume that a is a symmetric, continuous, and H- elliptic 
bilinear form. Let A be the operator associated with a and let oj, a be the 
constants from (2.2). Then 



(i) 
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where ch and cy are the operator norms of the embeddings V H and 
H V' , respectively. 

If the embedding of V in H is even compact, we write V H, then the 

c 

symmetric form a can be diagonalized in the following way. 

Theorem 2.2 (Spectral representation of symmetric forms). Let H and V 
be infinite dimensional Hilbert spaces with V H. Let o: V x V — > R be 

c 

a symmetric, continuous and H -elliptic bilinear form. Then there exists an 
orthonormal basis {e n : n G M} of H and A n G R mi/i Ai < A2 < . . . , 
limjj^oo A n = 00, suc/i that 

00 

V = \u G H : x n{u, e n ) 2 H < ooj, 

n=l 

00 

a(u, v) = ^2 ^(u, e n ) H (v, e n ) H (u, v G V), 

71=1 

00 

-4^ = X! K(u,e n ) H e n (ueV), 

n=l 

00 

£>(A) = [u G H : A «( u - e «)tf < 

n=l 

oo 

(u G 

n.=l 

This is a direct consequence of the spectral theorem for compact, self-adjoint 
operators (cf. [AU10, Satz 4.40 and 4.43]). The orthonormal basis in Theo- 
rem 2.2 is not unique, of course. But if u G D{A) is an eigenvector of A, i.e. 



4 



Au = An for some A G R, then u G span{e n : X n = A} = ker(A — A). Thus 

A n is in fact the n-th eigenvalue of A counting multiplicity. 

Next we extend our framework to non-autonomous problems. Let V, H 

be Hilbert spaces with V '—t H. 

d 

Definition 2.3. A Lipschitz continuous (non- autonomous) closed form is a 
function o: [0, r] x V x V — > R, where 

a) a(t, ■,■): V x V ->■ R is bilinear for all t G [0, r], 

b) there exists M > with 

|a(i, n, v) — a(s,n,n)| < M\t — s\\\u\\ v \\v\\ v (i, s€[0,t], B,1)6F), 

c) there exists M > with 

|a(t,it,u)| < M||tt|| v ||« || y (t G [0,r], u,t> G V"), 

d) there exist a > and to G R with 

Rea(t, u, u) + oj\\u\\ 2 h > a||«||y (t G [0, t],u G V). 

For the remaining section let a be a non-autonomous Lipschitz contin- 
uous closed form. For t G [0, r] we denote by A(t) the operator associated 
with the form a(t, •, •) and by A(t) the part of A(t) in H. We will need the 
following result from [ADL012]. 

Theorem 2.4. [ADL012, Theorem 5.1] Let B: [0,r] 6e strongly 

measurable such that 

(3\\uf H < Re(B(t)u,u) H < /3 _1 ||u||^ (t G [0,r], u e H) 

for a constant 1 > /3 > 0. Ze£ u° G V, / G L 2 (0,t; H). Then there exists a 
unique 

u G H l {<d,T;H)r\L 2 {{),T]V) 



satisfying 



B{t)u{t) + A(t)u(t) = f(t) a.e. 
u(0) = u°. 



Moreover, u G C([0,r];y) and 



1 1 2 1 1 • 1 1 2 

\ U \\l 2 (0,t,V) + IHIl 2 (0,t,H) - c 



where the constant C depends merely on f3,M,a,r and M. 
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The point of this theorem is that the solution u is in H 1 (Q,t;H) n 
L 2 (0,t;V); i.e. it has maximal regularity. This will be important in our 
context in order to make sure that the solution satisfies the desired boundary 
conditions. Moreover, it is remarkable that the solution is continuous with 
values in V. In our applications this implies in particular that the knot 
conditions are satisfied at each moment t > 0. 

Now we consider Theorem 2.2 in the non-autonomous case. Since the 
operators A(t) do not commute in general, the spectral decomposition (i.e. 
the orthonormal basis and the eigenvalues) will depend on t. We will show 
that the dependence of the eigenvalues is continuous in t. More precisely we 
assume the following. Let o: [0, r] x V x V — > K, be a Lipschitz-continuous 
closed form, where V H, and assume that a is symmetric. Consider the 

c 

operator A(t) on H associated with a(t, •, •) on H and denote by X n (t) the 
n-th eigenvalue of A(t) (counting multiplicity). 

Theorem 2.5. The function A n : [0, r] — > H is continuous for every n G IN. 

Proof. Let t G [0, r] and consider the constant oj of Definition 2.3 d). Then 
oj + A(t): D(A(t)) -)• H is bijective and (oj + A(t)) -1 G £(H) with (oj + 
A^^H = D(A(t)) C V. The mapping R: H -> V given by Ru = 
(oj+A(t))~ 1 u is bounded by Proposition 2.1. Since the embedding c. V — > H 
is compact by hypothesis, (oj + A^))^ 1 = to R(t) G C(H) is compact. 
Let s,t G [0,r]. Then 

||(o; + A(t))- 1 -(a; + A(5))- 1 || / . (fl) 

= \\(oj + A(t))~ l (oj + ^l(s) - (w + A(t)))(oj + ^(s))" 1 ^^) 
= || ( W + A(t))- 1 (A(s) - A(t))(oo + -4( S ))" 1 || £(H) 

< ||(^ + -A(*)) _:L llz:(v-',^)ll^ l ( s ) — ^-C*)!!^^,^')!!^ ^ •^ L ( s )) _1 ll J c(^ 5 v-) 

< M\t — s\ — , 

a a 

by Proposition 2.1 (ii), (iii) and property b) of the form a in Definition 2.3. 
Thus the mapping (oj + A(-))^ 1 : [0,r] — > C(H) is continuous. Now the 
theorem follows by the result in the Appendix (Theorem 6.1). □ 

Remark 2.6. Note that we could prove the continuity of the mapping (oj + 
J 4(-))~ 1 above without using property b) of the form a in Definition 2.3. We 
merely need that A(-) : [0, oo) —> C(V, V') is continuous instead of Lipschitz 
continuous. 

3 Diffusion in networks 

We apply the theory from the previous section to treat non-autonomous 
diffusion processes in a finite network. The network is represented by a 



6 



finite, simple, connected graph G with m edges ei , . . . , e m and n vertices 
vi, . . . ,v n . We identify all edges with the interval [0, 1]. We further assume 
that all the vertices have degree at least 2, i.e., that each vertex is incident 
to at least 2 edges. 

The structure of the network can be described by various graph matrices. 
Here we use the n x m incidence matrices $ + := (<frfj), &~ '■= an d 
<!> := (cf>ij), defined by 

6+ ■= I X ' if Vi = ^'W' 67. ■= I lj if Vi = e ^°^' 
* J ' | 0, otherwise, ij ' ' 1 0, otherwise, 

and 

Further, let T(vj) be the set of all the indices of the edges having an initial 
or endpoint at Vj, i.e., 

r ( v i) : = {j 6 {1, • • • , m} : <j>ij / 0} , i = 1, . . . , n. 

We will consider functions which are defined on the edges of the graph. 
These functions will be assumed to be continuous on the graph. By this we 
mean the following. Let / G C[0, l] m , f = (f 1 ,... f m ). For j G {1, . . . , m} 
we set /j(vj) := /,(0) if e 3 -(0) = v 4 and / 3 -(vj) := /j(l) if e 3 (l) = Vj. Then 
we say that / is continuous on the graph, if 

/j(vi) = /j(vi) (i^Gr(vi)) (3.1) 

for alH = 1, . . . ,n. 

Lemma 3.1. A function f £ C[0, l] m is continuous on the graph if and only 
if there exists a vector d G C™ swc/i i/iat ($~) T <i = /(0) and ($ + ) T d = /(l). 

Proof. Let / G C[0, l] m be a continuous function on the graph. For a vertex 
Vj denote by the "value of / at Vj", i.e. di := fjiyi) for any j G r(vj). 
Observe that for d := (<2j)i=i n we obtain (<l>~) T (i = /(0) and ($ + ) T d = 
/(I)- 

Conversely, if {^~) T d = /(0) and ($+) T d = /(l) for some d G C n , then 
for every j G r(vj) we have either d { = Sr.di = fj(0) or d { = Sf-di = fj(l), 
so (3.1) holds. □ 

We consider the following diffusion process in the network. At first we 
consider a finite time interval [0, r] where r > 0. 

• On every edge ej, j = l,...,m, evolution is governed by the heat 
equation: 

u j (t,x) = c j (t)-u'-(t,x) (t G [0,r], x G (0, 1)), (3.2) 

where c\(t), . . . , c m (t) are time-dependent positive diffusion coefficients. 
Here the dot refers to the derivative with respect to the time, and the 
prime refers to the derivative with respect to the space variable. 
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• We impose that for each t > the function u(t, •) = (ui(t, ■),..., u m (t, •)) 
is continuous on the graph; i.e., (3.1) is satisfied. 

• Boundary conditions at the vertices Vj, i = 1, . . . , n, are Kirchhoff-type 
conditions of the form: 

m 

£0ijMj(*K(*,Vi) =0 {t G [0, r]), (3.3) 

3=1 

where fJ>i(t), . . . , n m (t) are given time-dependent positive conductivity 
factors. Here we denote by Uj(i,Vj) the left derivative of the function 
itj(i, •) at if ej(0) = Vj and the right derivative at 1 if ej(l) = Vj. 

• The initial conditions on the edges ej, j = 1, . . . , m, are 

Uj(0,x) = u°j(x) (a? 6 (0,1)). (3.4) 

Before identifying the problem with an abstract Cauchy problem we 
write the continuity and boundary conditions in a more compact form. To 
this end we introduce weighted incidence matrices ^J(i) := UJ^(i)J and 

$~(t) := (uj~j(t) \ with entries 

+ ( U W, if V, = 6,(1), ^ ;= f K(t), if V, = 6,(0), 

U, otherwise, JJ U, otherwise. 

With these notations, the Kirchhoff law (3.3) can be written in matrix form 



as 



$-(t)-u'(t,0) = <f>+(t)-u'(t,l) (te[0,r}). (3.5) 

In order to use the previously introduced form methods we now first 
define the appropriate Hilbert spaces 

H := L 2 (0, l) m and V : = {/ G iT^O, l) m : / is continuous on the graph} 

equipped with the inner products 

rn .i 

(f,9)v = J2 (fj(x)g j (x) + f' :j (x)g' J (x))dx (f,g€V). (3.7) 

Recall that H 1 ^, 1) <^ C[0, 1] and so V is a closed subspace of ff^O, l) m . 
Note that V is a dense subspace of -£/\ For t G [0, r] we define the function 
a(V,-) : ^xF^Eby 

m r i 

a(t,f,g)=J2»j(t) ffaWjix) dx. (3.8) 

3=1 ° 
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Since (J,j(t) > for all j = 1, . . . ,m one easily gets that a(t, •, •) is a sym- 
metric, continuous, ^/-elliptic bilinear form. 

Proposition 3.2. Let t > 0. The operator associated with the form a(t, •, •) 
on L 2 (0,l) m is given by 

A(t) :=diagf-Mi(*)^2 J 

\ ' j=l,...,m 

D(A(t)) :={ueVn (H 2 (0, 1)) " : *~(t) • u'(0) = *+(t) • u'(l)} . 

Proof. Remember that the operator associated with the form a(t, •, •) 
on H is defined by 

D(A(t)) = {ueV :3he H s.t. a(t,u,«) = (h,v) H for all u G V}, 
A(t)u = /». 

We show that = A(t). 

First we prove A(t) C A(t). Let it G D(A(t)) and /i := ^4(t)u. Since 
u G V n (H 2 (0, l)) m , integration by parts yields 



a(t,u,v) = '*Tt*j(t) u 'j( x ) v 'j( 

7=1 1/0 



x) dx 



= ^^(t)[^(lH(l)-^(0H(0)j -$>j(*) / 4(x)^(x)dx 
i=i 3=1 

for all i> G V. By Lemma 3.1, for every v £ V there exists cZ„ G C n such that 
u(0) = ($~) T d„ and t> (1) = (& + ) T d v . Note also that the incidence matrices 
3> _ and 3> + have exactly one nonzero entry in every column (corresponding 
to the two ends of an edge) and have the same zero-pattern as the weighted 
incidence matrices 5>^(t) and 3>i(i), respectively. Thus 



- «i(0)«i(0)] = dj • (*+(*) • u'(l) - *~(t) • u'(0)) = 

3=1 

(3.9) 

for all v £ V. Hence a(t, u, v) = (h, v)h for all v G V. 

To show that A(t) C A(t), let u G and set h := A(t)u. Since 

a(t, u, v) = (h, v)h for all t? G V and £>(0, l) m C V we obtain by integration 
by parts as above that u G H 2 (0, l) m and hj = —fj,j(t)u'j. Thus 

m 

J2h(*) [^(iH(i)-^(oH(o)] =0 

3=1 
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for all v £ V. Using Lemma 3.1 as before we obtain (3.9) for every v G V . 
The entries of (d v )i are the joint values of Vj attained at vertex Vj for j G 
r(vj). Since this holds for arbitrary v G V, it follows 

*-- U / (0) = *+.«'(l). 

Hence u G D(A{t)) and A(t)u = h. □ 

Note that the domain D(A{t)) consists of continuous functions on the 
graph that satisfy the non-autonomous Kirchhoff-knot conditions (3.3). Thus 
in view of Proposition 3.2 we may rewrite (3.2)-(3.4) in the form 

tt(i) + A(t)u(t) = a.e. 
u(0) = u°, 

if fij(t) = Cj(t) for all j = 1, . . . , m. In order to treat the case \ij ^ Cj we 
introduce the operator 

B{t) := diag(6j(t)) G C(H) with 6j(t) := (3.10) 

Cj{t) 



Then (3.2)-(3.4) can be written as 

B(t)ii(t) + A(t)u(t) = (a.e. t G [0, r]) 
u(0) = u°. 



(3.11) 



Assuming some regularity properties on the time-depending coefficients c~(t) 
and fij(t) appearing in the equations (3.2) and the boundary conditions (3.3) 
we finally obtain the following well-posedness result of problem (3.2)-(3.4). 

Theorem 3.3. Let the functions fij : [0, r] — > R + be Lipschitz continuous 
and the functions Cj : [0, r] —> K, + be measurable. Suppose there exists an 
e G (0, 1) such that 

s < fJ,j(t), Cj (t) < e" 1 (J = l,...m, t G [0,r]). 

Let u° G V and F G L 2 (0, r; H) be given. Then there exists a unique function 

u G fl 1 (0,T;fl)nC([0 l T];7) 

siic/i i/ia^ /or all j = 1, . . . ,m the following holds 

Uj{t,x) = CjUj(t,x)+Fj(t,x) (t G [0,r], rr G (0,1)), 

m 

X! 4>ikVk(t)u' k {t,Vi) = (i G [0,r], i = 1, . . . ,n), 

fe=i 

nj(0, x) = u°Ax). 
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Proof. We check the conditions of Theorem 2.4. First we have to show 
that the function o: [0, r] x V x V — > R defined by (3.8) is a Lipschitz 
continuous closed form. This is true for M := max ig [ 0)T ] j=i.... ifn fJ>j(t), M := 
maxj = i v .. i?n Lj where Lj are the Lipschitz constants of the fuctions fij and 
u, a := e. 

The function B: [0, r] — > C(H) is measurable since the functions bj(-) : = 
are measurable. Moreover note that e 2 < < £~ 2 for all i G [0, r], 
thus we obtain that 

e 2 \\g\\ 2 H < Re(B(t)g,g) H < e~ 2 \\g\\ 2 H (3.12) 

for all g G H and t G [0, r]. Now given F G L 2 (0,t;H) we define F(t) : = 
B(t)F(t). Then by Theorem 2.4 we find a unique solution u of 

B(t)u(t) + A(t)u(t) = F(t) (a.e. i G [0, r]), n(0) = -u . 

Dividing by 6i(t) we see that u is a solution of the system as formulated in 
the Theorem. □ 

Thus we have well-posedness for every choice of the diffusion coefficients 
Cj(t) and the conductivity factors Hj(t) (up to slight regularity assumptions). 

Our aim is to examine stability. Therefore we extend Theorem 3.3 to the 
half line [0,oo). This is possible since r > was arbitrary. We formulate 
the result using an abstract (but equivalent) notation, which will be more 
convenient in Section 5. 

Corollary 3.4. Let the functions fij : [0, oo) —> R + be locally Lipschitz con- 
tinuous and the functions Cj : [0, oo) — > R + be measurable. Suppose there 
exists an e G (0, 1) such that 

e < fij(t), cj(i) < e^ 1 (j = 1, . . . m, t G [0, oo)). 

Then for every u° G V and F G Lf oc ([0, r); H) there exists a unique solution 
u G #/ oc ([0, oo); H) n Lf oc ([0, oo); V) n C([0, oo), V) of the non- autonomous 
Cauchy problem 

B(t)u(t) + A(t)u(t) = F(t) (a.e. t G [0, oo)) 
u(0) = u°; 

i.e. n|[ 0iT ] G L 2 (0, r; V) n -ff 1 (0, r; H) n C([0, r], V) is the unique solution of 

B(t)ii(t) + A(t)u(t) = F(t) (a.e. t G [0, r]), u(0) = u° 
for all t > 0. 
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4 Positivity 

In this section we show that the solution in Theorem 3.3 is positive whenever 
the initial value and the inhomogeneity F are positive. We start by a result 
which is of independent interest. 

Let (Q, E, /i) be a measure space. Then / G L 2 {Q) may be decomposed in 
its positive and negative parts f = f + — f~ where f + (x) = f(x) if f(x) > 
and f + (x) = otherwise, while / _ = (— /) + . 

Proposition 4.1. Let u G H l (0, r; L 2 (ft)). Define u+ : [0,r] ->• L 2 (0) &y 

u + (t):=u(t)+ (t€[0,r]). 
T/ien u+ G ^(0, r; L 2 (0)) and 

(u+)-(t) = u(t) • l { „(t)> } (o-e. t G [0, r]). (4.1) 
iJere l{ti(t)>o} z' s characteristic function of the set {x G f2 : > 0}. 

Proof. The mapping / H> / + : L 2 (f2) — > L 2 (0) is Lipschitz continuous. It 
follows from [AD013, Theorem 5.3] that u + G H l (0, r; L 2 (Q)). In partic- 
ular, it and ii + are differentiate outside a null set M C [0, r]. We assume 
that 0, t G M. Let t G [0, t]\M. We show that (4.1) holds for t G [0, t]\M. 
Let /i n 4 as n — > oo. Passing to a subsequence if necessary we find a null 
set N C f2 such that 

^(u(t ± /^(x) - u(t)(x)) -> u(i)(x) 

u(t ± h n )(x) — > u(t)(x) 
as n — t- oo for all a; G O \ iV. Now let x G \ N and consider three cases. 
1st case: u(t)(x) > 0. Then u(tzt h n )(x) > for n sufficiently large. Thus 
(u+)Xt)(x) = Jim g^(n+(£ ± h n ){x) - u + {t){x)) = u(t)(x). 

2nd case: u(t)(x) < 0. Then (« + )'(i)(x) = by the same argument. 
3rd case: u(t)(x) = 0. Then 

(u+)Xt)(x) = Jim ^(u+(t ± fc n )(x)) = 
since h n > and — /i n < 0. 

□ 
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Proposition 4.1 is somehow analogous to the classical result, valid for 
each open set Q C E d : If v G H 1 ^), then v + G H 1 ^) and 

DjV + = l {v>0} DjV (j = 1, . . . , d). 

Now we can prove the result on positive solutions. We consider the 
situation of Theorem 3.3. Let V and H be the Hilbert spaces defined in 
Section 3. For v G H we define 

v + : = («+)j=i,..., m 

and write u > if v = v + . Further we set v~ := (— t>) + , then v = v + — v~ 
and v + ,v~ > 0. Note that if v G V, then also "y + , f ~ £ V and 

«)'("7)' = ° 0" = 1, •••,"»). (4.2) 

Finally we denote the positive cones by 

H + := {u G H : u > 0} and V+ := {v G V : v > 0}. 

Theorem 4.2. Assume that all the conditions of Theorem 3.3 hold and ad- 
ditionally suppose that the functions b{ are Lipschitz continuous, i.e. bi G 
W 1,oo (0,oo). Moreover, let F G L 2 (0,t;H + ) and u° G V+. Then the solu- 
tion u provided by Theorem 3.3 is positive, i.e. it takes values in V+. 

Proof. Let c := ^ max i|l^lloo' wnere e > is as in Theorem 3.3. We set 
v(-) := e~ c 'u(-) and F(-) = e~ c 'F{-). Then F is still positive and, since u is 
a solution of (3.11) we have 

B{t)v{t) = e~ c \B(t)u(t) - cB{t)u{t)) 

= F(t) - A(t)v(t) - cB(t)v(t). 

Observe that v~ (0) = since n(0) > 0. Using (3.12) and Proposition 4.1 we 
obtain the following estimate for each t G [0, r]: 

e 2 \\v-(t)\\ 2 H < \\B^ 2 (t)v-(t)\\ 2 H = \\B 1/2 (t)v- (t)f H - WB^^v- (0)\\ 2 h 

= -2 [\b(s)v(s),v-(s)) h ds+ [\b(s)v-(s),v-(s)) h ds 
Jo JO 

= -2 [\f(s)-A(s)v(s),v-(s)) h ds 

JO 

+ f ((B(s)-2cB(s))v-(s),v-(s)) H ds 
Jo 

< -2 [ (F(s),v-(s)) H ds + 2 [ a(s,v(s),v-(s)) ds 

Jo Jo 

<2 / a(s, v(s), v~(s)) ds — 2 / a(s , v + (s) , v~ (s)) ds 
Jo Jo 

< -2 [ a(s,v~(s),v~(s)) ds < 0. 



Here we used that a(s,v + ,v ) = for all v G V by (4.2). This shows that 
v~{t) = e~ ct u~ (t) = 0, hence u(t) > for every t G [0, r]. □ 
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5 Stability 



Let V and H be the Hilbert spaces denned in Section 3. Note that V H 

c 

since ff 1 (0, 1) ^ L 2 (0, 1) (see [Brell]). Recall also the definition of the 

operator families (A(t)) t& m t00 \ and (-B(i))te[o,oo) from Section 3. We will 
assume all the conditions of Corollary 3.4. Note that then there exists a 
(3 £ (0, 1) such that 

P\\9\\n<\\B 1/2 (t)9\\H<0- 1 \\9\\ H (5-1) 

holds for for all g G H and all t £ [0, oo), cf. (3.12). Given some F £ 
Li oc ([0, oo); H) and u° G V we denote by u the solution to the abstract 
Cauchy problem 



(CP) 



B(t)u(t) + A(t)u(t) = F(t) (a.e. t £ [0, oo)) 
«(0) = it , 



given by Corollary 3.4. In this section we discuss the asymptotic behavior 
of u for several cases of the operator family (B(t)) t ^ oo y 

Since V *-¥ H we can apply Theorem 2.2 to the form a(t, •,•) for any 
t £ [0,oo). Thus the spectrum of the operator A(t) for any t £ [0, oo) 
consists of an increasing sequence of eigenvalues \i(t) < \2(t) < ••• with 
the corresponding orthonormal basis {e n (t) : n £ IN} of eigenvectors. The 
following properties of the first and second eigenvalue and the first eigenvec- 
tor of A(t) are essential for the analysis of the asymptotic behavior of the 
solution u. 

Lemma 5.1. Lett £ [0, oo). The normalized first eigenvector of A(t) equals 
e\(t) = -^7=1 and the first eigenvalue X\(t) = 0. The second eigenvalue \2(t) 
is positive. 

Proof. Note that (A(t)u, u)h = u, u) > for all t £ [0, oo) and all u £ 
D(A(t)), thus the first eigenvalue is non-negative by Theorem 2.2. Observe 
that ei(i) £ D(A(t)) and A(t)e\(t) = for all t £ [0, oo), which shows the 
assertions on e\(t) and \\(t). 

Since (e2{t) , e\(t)) h = we have that e^ (t) is not constant, i.e. ||e'2(i)||# 7^ 
0. Thus X 2 (t) = (A(t)e 2 {t),e 2 {t))H = a(t, e 2 (t), e 2 {t)) > by (3.8). □ 

Since the eigenvector e\(t) does not depend on t we set 

e x := -=1. 

Jm 



We can improve the result on pointwise positivity of A2(-)- 

Proposition 5.2. There exists a constant A 2 > such that X 2 (t) > A 2 for 
all t £ [0,oo). 
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Proof. By Theorem 2.5 the function A2 : [0, 00) — > [0, 00) is continuous. Sup- 
pose that there exists a sequence (i n )neiN C [0, 00) such that \2(t n ) — > 0. 
After taking a subsequence we may assume that t n — > t G [0, 00] as n — > 00. 
If t is finite X^ipri) ~ > ^2 (t) > as n — > 00 which is a contradiction. 

For the case t = 00 we also may assume that t n f 00 (after taking a sub- 
sequence). Recall the definition of the form in (3.8). Since the sequences 
(/Uj(i n )) ng ]N are bounded we may assume (after taking a subsequence) that 
they converge. Define continuous functions [O,^ 1 ] — > [e,e _1 ] with 
Mj(^n X ) = Mi(^n) where e > is chosen as in Corollary 3.4. For example 
define jlj affine on the intervals [in^n-i] an d set fij(0) = lim n _ 5 . 00 /j,j(t n ), 
then these functions are continuous on (O,^ 1 ] by definition and continuous 
in by a sandwich argument. For these functions we define the form a ana- 
loguos to a. We denote by A(t) the associated operator with the form a(t, ■, •) 
and by A(t) its part on H. The corresponding second eigenvalue of A(t) is 
denoted by \2(t). Now observe that ^(t" 1 ,?!, v) = a(t n ,u, v) and therefore 
^(^n 1 ) = ^2(t n ), hence ^(t^ 1 ) = \2(t n ) 0. On the other hand continuity 
of A2 (cf. Theorem 2.5 and Remark 2.6) implies ^(t^ 1 ) — > A2(0) > 0, which 
is a contradiction. □ 



5.1 Stability in the case B = Id 

In this subsection we consider the case that the operators B{t) in (CP) equal 
the identity operator for all t G [0, 00). Thus we have conservation of mass, 
i.e. (u(t),ei) H is constant. 

Proposition 5.3. Let u° G V , F G Lf oc ([0, 00); H) with (F(t), ei ) = for 
a.e. t > 0, and let u be the solution of (CP). Then 

(u(t), ei) H = (u°, e x ) H = -^(u , t) H (a.e. t G [0, 00)). 

Proof. Let t > 0. Then 

(u(t),ei) H - (v?,ei) H = I (u(s),ei) H ds = [ (F(s) - A(s)u(s), e\) H ds 

Jo Jo 

= f ((F(s), ei ) H -a(s,u(s),ei)) ds = 0. □ 
Jo 

We are now in the position to prove that in case B = Id the solution of 
(CP) consists of a constant part and exponentially stable part. 

Theorem 5.4. Let u° G V and F G L 2 (0, 00; H) with (F(i),ei) H = for 
a.e. t > 0. Then the solution u of (CP) can be decomposed as 

u(-) = u° + u(-) 

where u° := (u°, e\) H e\ = ^(u°, and u(-) := u(-) — u° is exponentially 

stable. 
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Proof. By Proposition 5.3 we have 

u(.) = u (.) - ( u °,ei) H ei = u(-) - (u(-),ei) H ei. 

Thus (u(-),e\) H = 0. This identity together with Theorem 2.2 shows that 

a(t,u(t),u(t)) > X 2 (t)\\u{t)\\ 2 H (a.e. t G [0,oo)). 

Let e > 0. Note that u is a solution of (CP) with initial value u° — vP instead 
of u°. Now by the product rule, Young's inequality and the above estimate 

l|fi(*)ll!r-||fi(0)l& 

= / 2(u(s),u(s)) H ds 
Jo 

= 2 [ (F(s)-A(s)u(s),u(s)) H ds 
Jo 

= 2 [ (F(s),u(s)) H ds-2/ a(s,u(s),u(s)) ds 
Jo Jo 

<h f\\F(s)\\ 2 H ds + 2e f\\u(s)\\ 2 H ds-2 f \ 2 (s)\\u(s)\\ 2 H ds 
Jo Jo Jo 

= i - 2 f\\ 2 (s) - e)\\u(s)\\ 2 H ds. 

Jo Jo 

By Gronwall's inequality we obtain 
\Wt)\\ H <(h -u% + i j)\F(s)\\ 2 H ds) 1/2 e W (- jT'cAaOO-e) ds) 
< (||n -n°||^ + ^ | o i ||F( S )||^d S ) 1/2 exp(-t(A 2 - e )) 
for all t > 0. This shows exponential stability of u if we choose e < A 2 . □ 

Remark 5.5. If the function F is only in Lf oc ([0, 00); ff) with ||i ? || i 2( t-H) — 
Ce^-V for some 6 > 0, C > and all i > 0, then n is still exponentially 
stable. 

Observe that the speed of the convergence of the exponentially stable 
part depends on the value of A 2 which is strongly related to the structure 
of the network (for a discussion on these connections we refer to [KMS07, 
Section 5]). 

5.2 Stability in the case bi < 

In this subsection we consider the case that the functions bi are Lipschitz 
continuous, i.e. bi £ W 1,oc (0, 00), and that bi(t) < for a.e. t S [0, 00). Again 
we denote by u the unique solution of (CP) . Conservation of mass does not 
hold in this situation, i.e. (u(-),e\) H may vary in time. Thus we cannot 
expect the same result as before but we can still show a nice asymptotic 
bahaviour of the solution. 
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Theorem 5.6. Let u° £ V and F £ L 2 (0,oo;H) with (F(t),ei) H = for 
a.e. t > 0. Then the solution u of (CP) can be decomposed as 

u = U\ + u, 

where u\(t) := (u(t),e\) H ei = ^(u(t),l.) H t converges to an equilibrium in 
H as t — > oo and u(-) := u(-) — u\(-) behaves exponentially stable. 

For the proof we need the following lemma. 

Lemma 5.7. Let f : [0, oo) — > [0, oo), g G L 1 (0,oo) such that for s < t 

f(t)-f(s)< f g(r) dr. 



Then the limit lim^oo f(t) exists. 

Proof. Define a n = inf{/(t) : t > n}. Then 

ai<02 <•••</(())+ / \g(r)\ dr, 

Jo 

hence a := limn^oo a n exists. Let t n > n such that f(t n ) < a n + -. Then 
for t > t n we have 

rt 1 />oo 

an < /(i) < /(*n) + / <7(r) dr < a n + - + / | 5 (r)| dr. 

Hence lim^oo f(t) = a. □ 

Proof of Theorem 5.6. Theorem 2.2 shows that 

a(t,u(t),u(t)) > X 2 (t)\\u(t)f H 

for a.e. t G [0, oo). Let e > 0. Now by (5.1), the product rule, Young's 
inequality and the above estimate 

p 2 \\ u (t)f H - r 2 \Ho)\\ 2 H < p 2 \Ht)f H - r 2 \Ho)\\ 2 H 

< \\B 1 / 2 (t)u(t)\\H-\\B 1/2 (0W n 



H 

= 2 / (B(r)u(r), u(r)) H dr + / (B(r)u(r), u(r)) H dr 
Jo Jo 

< 2 / (F(r) - A(r)u(r),u(r)) H dr 
Jo 

= 2/ (F(r), u(r)) H dr — 2 / a(r, u(r), u(r)) dr 
Jo Jo 

= 2 / (F(r), u(r)) H dr — 2 / a(r,u(r),u(r)) dr 
Jo Jo 

<h [ t \\F(r)f H dr + 2e f\\u{r)f H dr - 2 f X 2 (r)\\u(r)\\ 2 H dr 
Jo Jo Jo 

= h f\\nr)\\ 2 H ^ - 2 f\\ 2 {r) - s)\\u(r)f H dr. 



(5.2) 
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By GronwalFs inequality we obtain 



+ W J, W F ^\h dr ) ex P (t J, ~ £ ) dr 

«/ 



1/2 

exp^-^(A 2 -ejJ 

(5.3) 

for all t > 0. This shows exponential stability of « if we choose e < A 2 . 

Let < s < t. If we replace the lower integral limit in (5.2) by s we 
obtain 

WB^itHml- \\B^(s)u(s)\\h 

< Te /V«|Ih dr - 2 f(\ 2 (r) " e)||fi(r)||| dr. 

We insert (5.3) in the above estimate and conclude that ||i? 1 / 2 (t)tt(i)||^- 
converges for t — > oo by Lemma 5.7. Since the functions bi(t) are monotone 
and bounded they converge as t — > oo, thus also -B 1 / 2 (i) converges as t — > oo. 
Now 

(u(t),e 1 ) 2 H \\BV\t)e 1 \\ 2 H + \\BV 2 (t)m\H = 1 1 ^ 1/2 )^(* ) 1 1 S,r 

together with (5.3) shows that (u(t), e\)\ converges as t — > oo. Finally, since 
the function (u(-),ei) H is continuous also (u(t),ex) H converges as t — > oo, 
thus iii(t) = (u(t), ei) H e\ converges in H as t — > oo. □ 

5.3 Stability in the case fej(t) < 2cbi(t) for a c < A 2 

In this subsection we consider the case that the functions hi are Lipschitz 
continuous, i.e. 6j G I / t /1,oo (0, oo), and that 6j(i) < 2cbi(t) for some c < A 2 
and almost all t € [0, oo). Again we denote by u the unique solution of (CP). 
In this situation we can still decompose the solution as u = u\ + u and show 
that the function u is exponentially stable and u\ is a multiple of the first 
eigenfunction. But we do not know further properties of U\ like convergence 
as t — > oo, even boundedness is not clear. 

Theorem 5.8. Let u° £ V and F G L 2 (0,oo;H) with (F(t),ei) H = for 
a.e. t > 0. Then the solution u of (CP) can be decomposed as 

u = U\ + u, 

where u\(t) := ^(u(t),t) H t and u(-) is exponentially stable. 

Proof. We set v(-) := e~ c 'u(-) and F(-) = e~ c 'F(-), then since it is a solution 
of (CP) 

B(t)v(t) = e- ct (B(t)u(t) - cB(t)u(t)) 
= F(t) - A(t)v(t) - cB(t)v(t). 
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Let v(-) := v(-) — {v{-),e\) H e\. Theorem 2.2 shows that 

a(t,v(t),v(t))>Ht)\\m\\ 2 H 

for a.e. t S [0,oo). Let e > 0. Now by (5.1), the product rule and Young's 
inequality, we obtain for t > the estimate 

p 2 \\m\\ 2 H - r 2 Ho)f H < p 2 \nt)f H - 0-*\ko)\\h 



2 [\b(s)v(s),v(s)) h ds + [ (B(s)v(s),v(s)) H ds 
Jo Jo 

ft rt 



= 2 / (F(s)-A(s)v(s),v(s)) H ds+ / ((B(s)-2cB(s))v(s),v(s)) H ds 
Jo Jo 

<2 [ (F(s),v(s)) H ds-2 [ a(s,v(s),v(s)) ds 
Jo Jo 

= 2 [ (F(s),v(s)) H ds-2 [ a(s,v(s),v(s)) ds 
Jo Jo 

<h [ t \\Hs)\\ 2 H ds + 2e f t \\v(s)\\ 2 H ds-2 f X 2 (s)\\v(s)\\ 2 H ds 
Jo Jo Jo 

= i l*\\F(8)\\ 2 H ds - 2 f\x 2 (s) - e)\\v(s)\\ 2 H ds. 
Jo Jo 

By Gronwall's inequality we obtain 



(w\\ u °Wh + ^ f \\Hs)\\ 2 h dsf /2 exp f\x 2 (s)-e) ds) 
{w\\ u °Wh + W f\\n s )W 2 H d S ) V2 exp (- |(A 2 - 

J 



for alU > 0. Thus 

\\u{t)\\ H < (ji\\A\ 2 H + ^ f\\F( S )\\ 2 H d S ) V2 exp (- i(A 2 - c - 

J 

for all t > 0. This shows exponential stability of u if we choose e < X 2 — c. □ 



6 Appendix 

Let H be a separable complex Hilbert space of infinite dimension. Suppose 
A is an unbounded self-adjoint operator on H with a (A) C [coo, oo) and 
compact resolvent (i.e. R(u,A) = (u — A)~ l is a compact operator on H 
for some uj < ojq). Then by the Spectral Theorem, the space H has an 
orthonormal basis {e^ : k G IN} such that e& G -D(-<4)> = A^e/j where 
wo < Ai < • • • < A,t < Afe + i < . . . and lim^oo A& = oo. As a consequence 
{Afc : k £ IN} is exactly the set of all eigenvalues of A and hence A& is the 
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k-th eigenvalue of A counting multiplicity. For n £ M let A n be a further 
unbounded self-adjoint operator on H with a(A n ) C [ujq,oo) and compact 
resolvent. We denote the k-th eigenvalue of the operator A n by AjJ. 

Theorem 6.1. Suppose \\R(u), A n ) — R(u, A)\\£( H \ — > as n — > oo for some 
to < ojq. Then lim. n -+ rx A£ = A& for every k £ IN. 

We show two auxiliary results before proving the theorem. 

Lemma 6.2. Let loq < /ii < • • • < a p < u p+ i < . . . be all different eigen- 
values of A, Hp having multiplicity m p £ IN, and let ujq < < A^ +1 < . . . 
be such that the following two conditions hold. 

a) For each p £ I, e > such that \i p -\ + £ < fi p < Hp+i — £ there exists 
a no £ IN such that for all n > no one has 

#{k : Afc £ [fi p -£,H P + e]} = m p . 

b) For all k £ IN each limit point of (A^) n6 iN is in {u p : p £ IN}. 

Then lim n _ s>00 A^ = Xk for all k £ IN where Xk = /i p /or M p _i < k < M p , 
M p = m\ + • • • + m p for p £ IN and M$ = 0. 

Proof. The proof is given by induction over p. Let p = 1 and let e > such 
that Hi + e < fj,2- By a) there exists a uq such that A^ fi < /^i + £. Thus for 
k = 1, . . . , Mi, the sequence (A^) n6 ]N is bounded and by b) it has u\ as the 
only possible limit point. Hence lim n _ > . 00 

Xl = ni for = 1, . . . , Mi. 

Now let p £ IN, p > 2 such that the assertion of the lemma holds for 
p — 1. Let e > such that Hp-i < /%> ~~ £ < /%> < + e < Mp+i- By 
the inductive hypothesis there exists no £ IN such that for all n > uq, 

_i — ^ A*p-i < Hp — £■ It follows from the assumption a) that there exists 
n\ > no such that X\ [p < n p -\-e for all n > n\. Hence the only possible limit 
point of (A£) ne iN is u p whenever M p _i < k < M p . Thus liuin^oc A^ = u p 
for M p _i < < M p . □ 

Lemma 6.3. Lei w < i^o such that ||i?(u;,^4 n ) — R{uo, A)\\^^ — > as 
n oo. If a £ /?(j4), i/ien i/iere exists a no £ IN such that /i £ p{A n ) for all 
n > no- 

Proof. Suppose that /i £ o~(A n ) for infinitely many n £ IN. After taking a 
subsequence we may assume that [i £ o~(A n ) for all n £ IN. Thus there exist 
eigenvectors e n such that A n e n = ne n for all n £ IN. We get 

e n = R(u,A n )(uj - A n )e n = R(u, A n )(u - n)e n 

= (uj - u)(R(u,A n ) - R(u,A))e n + (a; - h)R{lo, A)e n 

for all n £ IN. Since R(ui,A) is compact and ||e n ||^ = 1 after taking a sub- 
sequence we may assume that R(u,A)e n converges. This, the convergence 
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of the resolvent and (6.1) imply that also e n converges to some e as n — > oo 
and that e = (uj — fj,)R(uj, A)e or Ae = /xe. This contradicts the assumption 
that [i G p(A). We get that \i G cr(A n ) for only finitely many n £ IN, which 
proves the lemma. □ 

Proof of Theorem 6.1. We apply Lemma 6.2 to X%, /x p , m p defined above. 
Condition b) is satisfied by Lemma 6.3. Next we want to prove condition 
a). Let p G IN and e > such that p p ~i + e < p p < [i p +i — e. It follows 
from b) that there exists no G IN such that [i p -\ + e, p p +i — e G p(A) for all 
n > uq. For n > uq we denote by 

P n := J R(ui,A n ) du (6.2) 

the spectral projection associated with (/i p _i + e,u p +i — e) n a {A). Here V 
is the circle with center ^(u p _i + fi p +\) and radius r = ^(p p+ i — n p -i) — e. 
Similarly, let 

P := J R(u,A) du. (6.3) 

Then \\P n — P — > as n — > oo. Hence there exists n\ > no such that 
\\P n — P \\£fjj\ < 1 for all n > n±. Consequently P n and P have the same 
dimension [Kat66, I. Theorem 6.32]. This shows that condition a) holds. 
Now the assertion of Theorem 6.1 follows from Lemma 6.2. □ 
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